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Two models for the preheat stage of conventional liquid fuel droplets and of emulsified fuel droplets in 
combustion gases are analyzed theoretically. These models contain the effects of transient heat conduction to the 
droplets. In the first model, the droplet and gas temperatures var\ temporally hut only the gas temperature varies 
spatially; i.e., the droplet temperature is spatially uniform. Numerical examples, computed from this model, for 
both the droplet and gas temperatures are given. In the second model, both the droplet and gas temperatures vary 
spatially and temporally. Numerical examples computed from this second model for the surface and average 
temperature of the droplet are given. These analyses show that the temperature gradients inside droplets of oil 
and water are small compared to those in the combustion gases near the droplet; that temperature profiles given 
by both models are very similar. In particular, the predicted times at which micro-explosions are expected to 
occur agree within 10 percent of each other. 
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1. Introduction 

Before liquid fuels can burn effectively in combustion chambers, they must be reduced to minute 
droplets. This atomization of Liquid fuels may involve one or more of the following procedures: (1) mixing the 
fuel with high pressure air or gas; (2) injecting the fuel through small orifices at high pressures; and (3) more 
recently, superheating emulsified liquid fuels to their "micro-explosive" limit. Effective atomization occurs 
as follows. First, thin sheets of liquid fuel form. They then become unstable and form ligaments and large 
drops. Finally, the segments of ligaments and large drops disintegrate further to small droplets [I]. 1 

The ignition of the resulting droplet of conventional fuel consists essentially of two stages [2]. During the 
first preheat stage, heat flows from the hot surrounding gas to the droplet without significant evaporation. In 
the second stage, ignition occurs in the gaseous mixture surrounding the droplet and consisting of hot air 
(oxidizing atmosphere) and fuel vapors. The preheat stage is dominated by transient heat flow and is the 
subject of this paper. 

Several researchers have proposed that emulsified fuels may yield increased combustion efficiencies; 
decreased pollutant emissions, such as soot, CO, and oxides of nitrogen; and lower maintenance costs [3, 4]. 
Emulsified fuels consist of two or more components. At least one component has a much higher vapor 
characteristic than the primary fuel component. The component with the higher vapor pressure, for example 
water, is dispersed throughout the primary fuel. Such emulsified fuels may be generated by ultrasonic 
techniques. 

The corresponding preheat stage for emulsified fuels consists of two phases. For illustrative purposes, we 
consider in this paper a two component emulsified fuel. Injection of such an emulsified fuel into a combustor 
creates compound drops in which the primary fuel surrounds a number of smaller droplets of the dispersed 
component. Because the smaller droplets have a lower boiling point than that of the surrounding primary 
fuel, vaporization of the former occurs at superheat temperatures above the normal boiling point. Superheated 
boiling or spontaneous nucleation of the dispersed component is explosive and occurs at surprisingly 
reproducible superheat temperatures for ultraclean emulsions [5, 6]. The resulting sudden change in the 
volume occupied by the dispersed phase produces catastrophic shattering of the primary fuel drop into a 
number of much smaller droplets (fragments). This then is the first phase of preheating. The second phase of 
preheating for emulsified fuels occurs when additional heat transfers to these smaller droplet fragments 
without significant evaporation and before ignition occurs. 



1 Figures in brackets indicate the lite* rature references at the end of this paper. 
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In this paper, we develop two models for predicting the transient heat flow to conventional fuel droplets 
undergoing preheat and to emulsified droplets undergoing the first phase of preheat. We shall use these 
models to compute the dependence of the temperature of the droplet and its surrounding, hot oxidizing gas 
(e.g., air) upon distance from the droplet center and time from injection into the hot oxidizing gas. 

In the first model, the droplet and gas temperatures vary temporally but only the gas temperature varies 
spatially. That is, the droplet temperature is spatially uniform. We shall refer to this model as the zero- 
gradient model. Some researchers propose that in many situations the heat transfer rate within the droplet is 
much faster than is possible with conductive transfer alone and that internal circulation may be sufficient to 
maintain a spatially uniform temperature within the droplet [7]. Hence, the zero-gradient model of this paper 
is an appropriate one for such cases. In the second model, both the droplet and gas temperatures vary 
spatially and temporally. We shall refer to this second model as the finite-gradient model. 

We include these two models for the following reasons. First, the numerical procedures to evaluate the 
temperature profiles predicted by the zero-gradient model require an order of magnitude less computer time 
than the numerical procedures to evaluate the temperature profiles predicted by the finite-gradient model 
require. We expect that the zero-gradient model would be a particularly useful one for any future theoretical 
analysis of the ignition stage of the droplet. Second, the zero-gradient model readily gives the spatial 
dependence of the gas temperature near the droplet. We anticipate that such spatial details will be required 
if one were to consider dense sprays of individual burning particles. And third, we expect that the finite- 
gradient model is the more appropriate one for very short times prior to the micro-explosions of emulsified 
droplets and for those situations in which the thermal conductivity of the gas is comparable to that of the 
droplet. 

In the next section, we describe the physical assumptions for the preheat stage which are common to both 
models. Section 3 contains the mathematical description of these assumptions. We derive in section 4 
profiles of the droplet and gas which are predicted by the zero-gradient model. In section 5, we give the 
expressions for the surface temperature of the droplet and for the spatially averaged temperature within the 
droplet for the finite gradient model. We present in the last section numerical examples for oil-water 
droplets. We find that the reduced temperatures given by the two models are within 15 percent of each other 
for small reduced times and are within 5 percent of each other for large reduced times. 

2. Preheat Stage 

In models for droplet preheating, many researchers consider the droplet to be a sphere with temperature 
independent thermal properties and neglect diffusion effects [2]. Others include thermal diffusion effects but 
later in their theoretical analyses neglect some terms in the diffusion equations [8]. In addition, most 
theoretical investigators neglect the influence of gravity, forced convection, and radiant heat transfer. 
Experimental techniques such as freely falling combustion chambers can be used to eliminate gravitational 
effects and forced convection. Proper selection of experimental parameters such as ambient pressure and 
initial droplet size may minimize the other effects not included in the model for preheating. Several 
researchers have shown that except for heavy fuel oils, radiant heat transfer from the hot gas or from 
adjacent droplets is negligible [7, 9]. 

The models examined in this paper contain thermal diffusion effects. We assume that a spherical liquid 
fuel droplet of radius r d is inserted at time t = into a hot atmosphere (e.g., air). The initial temperature of 
the droplet at time t = is T^. The hot oxidizing gas is unbounded and at a constant temperature T g0 
infinitely far from the droplet. During the preheat stage, we assume that the radius of the droplet does not 
change. That is, no evaporation occurs. Evaporation will be considered in a subsequent paper. 

A conventional fuel droplet is essentially homogeneous and has well defined thermal properties. However, 
an emulsified fuel droplet composed for example of water and oil, is not microscopically homogeneous. In 
order to make the models for emulsified droplets amenable to calculation, we must make an additional 
assertion. For the purposes of transient heat flow analysis, we assert that representative values of the thermal 
properties of an emulsified fuel droplet exist and that they characterize properly the overall thermal behavior 
of the droplet. These representative values are essentially uniform throughout the droplet and are most likely 
bounded by the thermal properties of the fluid components in the emulsified droplet. 
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We denote the density, specific heat at constant volume, and thermal conductivity of the droplet by d d , 
C d , and k d , respectively and of the hot gas by d g , C g , and K g , respectively. We assume that these thermal 
properties remain spatially and temporally constant during preheat. In addition, the ambient pressure of the 
host gas, p g , does not change. 

Using the above models we seek to compute the dependence of the temperature T(r, t) upon the radial 
distance r and time t. In particular, we shall compute the time required for a conventional fuel droplet to 
reach a given temperature and for the emulsified fuel droplet to reach its superheat limit T SL ; i.e., the 
temperature near which the micro-explosion occurs. Because the adsorption of impurities at the interface 
between liquids might alter surface and interfacial tension, most water droplets in oil vaporize explosively 
somewhere between 240 and 270 °C [3, 5J. 



3. Theoretical Analysis 

We find it convenient at this point to introduce several dimensionless quantities. We denote the reduced 
temperature r9(Tj, t) by the expression 

0(i7, t) = [7V, t) - T M )/(T g0 - T d0 )l (1) 

The dimensionless distance is 7) = (r/r d ) and the dimensionless time is T = (t/t d ). In this paper, we shall 
denote dimensionless quantities by Greek letters. The time t d is proportional to the time required for a heat 
pulse to diffuse across the droplet and is given by t d — r d 2 /a d 2 where a d 2 = {K d /d d C d ) denotes the thermal 
diffusivity of the droplet. Even though the reduced time T may be mathematically unbounded, it does have a 
physically determined upper limit for the preheat stage. The physically meaningful values of T must be less 
than that value of T = T max for which 0(T/, T max ) corresponds to the boiling point of the liquid in the droplet. 
The Fourier heat conduction equations for the droplet and its surrounding £as describe the time and space 
dependence of the temperature when local thermodynamic equilibrium exists. For the case in which the 
isothermal surfaces are concentric spheres centered about the droplet center at r\ = 0, these equations are, 
forr > and 17 < 1 (inside the droplet) 

d 2 d 2 d6 d dO d 
drf j] or) or 

and for T > and 7] > 1 (outside the droplet) 

d 2 #g 2d0y = o2^, (3) 

drj 2 7) 3t) dr 

where 7) > and a = a d ja g . The quantity a g z = (K g /d g C g ) is the thermal diffusivity of the gas and the 
subscripts d and g refer respectively to the droplet and the gas. 

These Fourier heat conduction equations require a statement of the boundary conditions before solutions 
are uniquely defined. The reduced temperature 8 has the form 



and 



and it becomes for T > 



d (7], r) = for 7) < 1 and r < 



Og(v, t) = 1 for 7) > 1 and r < 






r, r) for 7) < 1 

*i.*M*/- ^ for , >L 
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The boundary conditions do not depend explicitly upon time. The temperature is finite everywhere 
namely, #d(0, r) is finite and 

Lim 6 9 (r), t) = 1. (4) 

Also, the condition that no heat flux exists at the center of the droplet becomes 

^L-OL (5) 

07) 

The continuity of the temperature across the droplet-gas interface is 

e a (h t) = e g (i, T ). (6) 

The conservation of heat flow at the interface gives, 

afldfr,T) | n dO g (v,r) \ 

|t,«i — p Jn-i> \<) 

07) 07) 

where /3 = K g /Ka is the ratio of thermal conductivities for the gas and droplet. 

Equations (2) through (7) represent the mathematical description of the above model for transient heat flow 
during the preheat stage of an individual fuel droplet injected at time t = into hot oxidizing gases. We 
solve the equations by taking their Laplace transforms with respect to time. The Laplace transform of the 
reduced temperature (tj, r) is denoted by 

c 

4>(i7, cr) = I exp (—err) 6(7), r) dr, 
Jo 

where <x is the dimensionless Laplace transform variable. The Bromwich integral 

t r€+i8 

6(7), r) = Lim ; exp (ar) 0(17, cr) da, (8) 

s-*oo 2m J € -j$ 

expresses the temperature in terms of the Laplace transform. The quantity € is chosen sufficiently large so 
that the integral, 



exp ( — ot)|0 (17, r)\dr 
Jo 



exists. 



Introducing the dimensionless variable X = cr 112 17, we express the Laplace transforms of eqs (2) through 
(7), respectively in the forms 

a 2 2 a 1 

$X X $X J 



— 2 + -— -« 2 K,(x) = — 



(10) 
J °" 

— -1 =0, (11) 

d X l x =o 
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im <t> 9 (x) = o- 1 (12) 

I)-»00 

<J>„( X = p) = <t> g (x = P) (13) 

^| =/A , (.4) 

a 'x=p d x >X=P 

where p = cr 1/2 . 

The solution to eq (9) and the condition (11) require that the Laplace transform for the droplet has the form 

<& d (x) = r d x- 1 sinh(x). (15) 

Similarly, with Xo = °0C^ tne solution to eq (10) and the condition (12) require that 

QgiXn) = JoXs,' 1 ^ X P (Xg) + P' 2 ' (16) 

The remaining continuity condition (13) and the conservation of heat flow eq (14) yield two simultaneous, 
inhomogeneous equations for the coefficients y^ and y g . Solving these equations, we obtain the respective 
Laplace transforms for the droplet and the gas; 

0> d (x) = A-^fl + ap) sinl. (x) for t, < 1 (17) 

and 

*,(x) = -A" 1 {p cosh (p) - sinh (p)} exp {- ap (t, - 1)} + p" 2 , (18) 

for 7] > 1, and where 

A = P 2 tj [{#1 + ap) ~ 1} sinh (p) + p cosh (p)]. (19) 

Inserting eqs (17) and (18) respectively into the Bromwich integral (8) gives us the reduced temperatures 
0<fO?» T ) f° r tne droplet and g {r), r) for the gas. In the following section 4, we give expressions predicted by 
the zero gradient model for the reduced temperature 6 d in the droplet and for the reduced temperature 6 g in 
the gas near the droplet. In section 5, we present the exact evaluation of the reduced temperature 0j inside 
and on the surface of the droplet for the finite gradient model. 

4. Zero-Gradient Model 

Two different but equivalent ways exist for obtaining the expressions which correspond to the zero- 
gradient model. We shall give both. 

First, for sufficiently large reduced times r, the major contributions to the Bromwich integral (8) arise 
from the region corresponding to values of |cr| = |p 2 | which are less than one. This suggests that we consider 
expansions in powers of p. We divide eq (17) and (19) by sinh (p) and for small values of p approximate p 
cosh (p) by 1 + (p 2 /3) and sinh (p) by p. For even smaller values of X = VP because t) < 1, we 
approximate sinh (x) by \ an d nn d that, 

•«fo) -♦ 0*(X) = 3/9(1 + ap)[p 2 (p + P + )(P + P-)T\ (20) 

where 

p ± = (3a / 8/2)[l ± {1 - (4/3/W 2 ]. (21) 

We apply the partial fraction method to eq (24) and rewrite it in the form, 

•*(X) = P~ 2 " [<P + P+)(P + P-K 1 - ( 22 ) 
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Observe that for sufficiently long reduced times T, the droplet reduced temperature da becomes independent 
of rj; i.e., it is spatially uniform and depends only upon the reduced time T. 
Procedures similar to that by which we obtain eq (20) give us that 

*«(X) -> ®,z(x) = P" 2 - blip + P+)(p + P-)V exp {-api-n - 1)}. (23) 

Using the relation that 

[(p + p + )(p + p_)]-> = (p_ - p + )-' {(p + p + )"' - (p + p.)-'} (24) 

and referring to page 1023 and page 1026 of Abramowitz and Stegun [10], for 7) < 1, we evaluate the 
Bromwich integrals for eqs (22) and (23). The results are, for 7) < 1, 

Odz(r) = 1 + (p+ - p_) _1 {-p+ exp {p + 2 T) erfc (p+ Vr) 

+ p_ exp (pJr) erfc (p_ Vr)} (25) 

and for 7] > 1 , 

Ogziv, r) = 1 + {7j(p + - P-)} -1 [-p+ exp{p+ 2 r + q^).(t7 - 1)} 

erfc { P+Vt + (0(19 - l)/2 Vr)} + p_ exp {pJr + op (r) - 1)} 

erfc{p_ Vr + (a(rj - l)/2\/r)], (26) 

where the subscripts <iz and gz denote the droplet and gas reduced temperatures for the zero-gradient model 
and where the complementary error function of complex argument z is denoted by erfc (z). Equations (25) 
and (26) are valid only when p + =^ p_ . 

For fuel droplets in combustion gases, the values of a and /3 are such that p+ equals the complex 
conjugate of p_; i.e., p + =p_*. We then analytically continue the functions appearing in eqs (25) and (26) 
to their values for complex arguments in order to compute the reduced temperatures. This is a tedious 
process and we refer the reader to [11] for the details. 

In order to compute the reduced times T for which eqs (25) and (26) correspond to valid inverse transforms 
of eqs (17) and (18), we consider the previous inequality cr 1/2 < 1. We divide both sides of this inequality by 
O" and then take the inverse Laplace transform. This yields the inequality (ttt)~ 112 < 1 or t > 77" -1 . 
Remembering that T = {t/ta), we have that for times t > (ra 2 /7rcLa 2 ) the temperature gradients in the droplet 
become negligible. The time (ra 2 /TTaa 2 ) corresponds to the time which is required for a heat pulse to diffuse 
a distance r^ in the droplet. 

There is a second way in which to obtain the same reduced temperatures given by eqs (25) and (26) 
without making any assertions about the value of T. The second argument involves considering the 
conservation of heat flow at the interface for the case in which the total thermal conductivity of the droplet 
due to both Fourier heat conduction and internal circulation exceeds that for the gas. As we mentioned in 
section 1, the internal circulation inside a droplet may be sufficient in many cases to maintain a spatially 
uniform temperature within the droplet. When such conditions prevail, the Fourier heat conduction equation 
for the droplet, eq (2), and the boundary condition (7), are replaced by a single heat transfer equation. For 
the case in which the droplet temperature is spatially uniform, Ta(r, t) = Ta(t) = T a (ra, t), the heat flux into 
the droplet, 

{(4TTr d 3 /3)ddC d (dT d /dt)/4TTr d 2 } (27) 

must equal the heat flux from the gas, K (dT g (r, t)ldr) r=rd . This last equality replaces eqs (2) and (7) and 
becomes in terms of the dimensionless quantities used in this paper, 



1 d6 dv „d 
3 rfT df) 
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(28) 



where 1 6 dv and 6 gu are the reduced temperatures for the spatially uniform case discussed above. The Laplace 
transform of eq (28) is 



**.(o-)= P — OgJiri, <r) 
3 dr) 



(29) 



and this eq (29) replaces eqs (9) and (14). Remembering that 4> dM (cr) = ^ u (l, O"), we use eqs (10) and (12) 
with <t> g replaced by <t> gu and eq (29) to obtain the result that <t> gu (j), cr) = ^ gz (x)^ where <f> 0Z ix) is given by 
eq (23). Hence, we then have that for all values of reduced time T ^ and when (d7^(r, t)/dr) = for all r 
< r , the reduced droplet temperature is #d U (r) = ^dz( T ) an d OguiVi T ) = OgziVi r )> where 6<i z (t) ar, d gz (r], 
t) are given by eqs (25) and (26). 

For the case of an oil or water droplet in hot air at atmospheric pressure, the quantity (4/3/3a 2 ) is of the 
order of 10 4 and the quantities p+ and p_ are approximated very well by p ± ~ ± t(3/3) ,/2 . Using these 
values in eqs. (22) and (23) and assuming that a(t) — 1) <^ 1 in eq (23), we obtain the rather simple 
approximations to eq (25) and to eq (26); namely, 

UmO dz (r) = 0«f O (r), 

ot->0 



and 



LimO gg (r), r) = gO (v^)^ 

ot-»0 



where 

M (r) = 1 -exp(-30r), (30) 

and 

e B0 {r,,T)= 1 -T,- 1 e Xp (-3/3r). (31) 

From eq (30), we observe that (3/3) -1 is the reduced time required for (1 — 840} to reach (e _1 ) times its value 
at r = 0. We shall see in section 6 that for oil and water droplets in air, the numerical predictions of eqs 
(30) and (31) agree very well with the numerical predictions of the exact equations (25) and (26). 

5. Finite Gradient Model 

In the finite gradient model, the reduced temperature distributions of the droplet and the surrounding gas 
are obtained by directly evaluating the Bromwich integrals, 

n re+id 

0(7), r) = Lim — ; exp (o-t)4>(t}, <t) da. (32) 

5-*» 27Ti ./e-iS 

<I>(T7, O"), or equivalently 4>(x) with \ ~ v V%r ), 1S given by eq (17) for the droplet and eq (18) for the 
surrounding gas. 

To evaluate the Bromwich integral it is first necessary to investigate the nature of the integrand exp 
(OTYf>(r), cr). No distinction will be made in this discussion between the droplet and the gas as the 
integrands in both cases exhibit the same general properties. Since exp (crr)<I>(i7, cr) is not an even function 
of X ^ 1S not a single valued function of (J, and it possesses a branch-point at cr = 0. The quantity exp 
(cr^)4>(T7, cr) also has a simple first-order pole with residue unity at cr = and singularities at the nonzero 
roots of A given in eq (19). Simple considerations of A indicate that no pure real or pure imaginary roots 
exist; but they do not rule out the existence of complex roots. We present in appendix A a general proof that 
the function A does not have any complex roots [12]. 
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We now evaluate the Bromwich integral. We compute the contour of equation (32) by two quarter circular 
arcs in the upper and lower left hand plane, by two parallel lines running above and below the negative real 
cr-axis, and by a small circular arc traversing a path around the origin cr = 0. The reduced temperature 
distributions for the droplet and the gas are then obtained by applying Caueky's Theorem and by letting 8 — ► 
00 in eq (8) and the radius of the smaller circle about the origin tend to zero. The quantity € is taken to be an 
arbitrarily small positive number. The integrations along the quarter-circular arcs vanish as their radii go to 
infinity. The simple first order pole at cr = contributes unity, and the two infinite integrals above and 
below the negative real cr-axis are evaluated with cr set equal to £; 2 i l7T and ^ 2 e~ l7r , respectively. We then 
obtain the following expressions for the reduced temperature distributions: 

n ( \ _ ■. . ?^ f" sin fr}(£ cos £ - sin g exp (-TJ?) dj 

W»' T > " + ^ J Jfp sin 2 f+(fcoB f- (1 _0 sin f )2 W 

and 



2_ f N(r,, fltf cos | - sin g exp (-rftdt 
irq Jo c?f?? sin 2 £ + {£ cos £ - (1 - j3) sin £) 2 



e B (r,,T) = i + - '^?\:r\ 7 ;r ' . ;^ (») 



where iV(i7, f ) = £ 1 sin (a(j) — 1)£)[£ cos £ — (1 — (3) sin £] + a/3 sin f cos (a(r) — l)f). The average 
temperature of the droplet is given by 



d (avg, t) = r) 2 6 d (r),T)dri/ if dr) 
Jo Jo 



6t>!j8 p r 2 (^ cos £ - sin £) 2 exp (- t£ 2 )<^ 
_1 77-1 a?/?f 2 sin 2 f + (f cosf -(1 - )8)sinf) 2 ' 



(35) 



and a similar expression for the surrounding gas may easily be obtained from eq (34). 

The above rather complicated integrals must be evaluated numerically. A short outline of the procedure 
which we used is given in appendix B. 

We discuss in the next section the representative values of the reduced surface temperature of the droplet 
6 8 (t) ■ g (l, r) given by eq (33) with 7] = 1 and the spatially averaged reduced temperature throughout the 
droplet flavW — 0<j(avg, r) given by eq (35). 

6. Numerical examples and Conclusions 

In this section, we give some illustrative numerical examples for the predictions made by the zero-gradient 
and finite gradient models. The input data for these calculations are the thermal conductivity, the density, 
and the specific heat of the fuel droplet and the combustion gas. Values for the above properties of a 
conventional fuel droplet are readily available. Table 1 contains typical values of medium weight fuel oils. 



TABLE 1. Thermal properties 

The thermal conductivity, the density, and the specific heat at constant volume are denoted respectively by K, d, and C. 

K d C 

(J/cm s K) (g/cm 3 ) (J/g K) 

oil 1.45 X 10" 3 0.95 2.24 

water 6.75 X 10 3 1.00 4.20 

eas 2.70 X 10" 4 1.30 X 10^ 3 1.00 



Values for the above properties of an emulsified fuel require more consideration. We expect that an 
emulsified droplet containing oil and water will have thermal properties which are bounded by those for the 
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pure fuel oil and pure water components. Hence, we list also in table 1 the values for the thermal properties 
of water. An emulsified droplet of 10 per cent water and 90 per cent fuel oil probably has an overall thermal 
behavior which is very close to that of fuel oil. Nonetheless, we include calculations for the pure water 
droplets to give bounds on the predictions and to test the sensitivity of the predictions upon the thermal 
properties of the droplet. Past experience indicates that whenever the thermal conductivity of the spherical 
droplet being heated exceeds the thermal conductivity of the medium surrounding the droplet, the thermal 
response of the droplet-gas system is dominated substantially by the thermal properties of the surrounding 
medium. We use the values of air at one atmosphere of pressure to represent the hot oxidizing combustion 
gases. These values are also given in table 1. Table 2 contains values for the two quantities a and /3. The 
radius of the droplet r<j is 50 /Jim for all the numerical examples considered in this paper. 

Table 2. Reduced parameters 

The dimensionless quantity <* 2 is the thermal diffusivity of the droplet #d = (KdJdaCd) divided by the thermal diffusivity of the gas 
a\ = {Kg/djCg) and the dimensionless quantity (3 is the thermal conductivity of the gas K g divided by the thermal conductivity of the 
droplet. The density and specific heat of the droplet and the gas are denoted respectively by d& and C^ and by d and C . 

a (3 

oil droplet 0.0572 0.186 

water droplet 0.0880 0.040 



Figures 1 and 2 compare respectively the numerical predictions of the zero gradient and finite gradient 
models for oil droplets in air. Figures 3 and 4 compare respectively the numerical predicitions of the zero 
gradient and finite gradient models for water droplets in air. From these figures, we observe that the four 
reduced temperatures 0d*( T )? 0do(T), Qs{ T )i an d #av( T ) all agree to within 15 per cent of each other for values 
of T ^ 0.5 and to within 5 per cent of each other for values of T « 2.5. Because all four quantities approach 
unity as T approaches infinity, the agreement among the values improves for values of r > 2.5. We also find 
that the water droplets require longer reduced times T to achieve reduced temperatures comparable to those 
of the oil droplet. 

Because the computer time required to evaluate 8 (t) and a v( T ) exceeds the computer time to evaluate 
0(1 2 (t) by at least a factor of ten and because the numerical agreement among 0^(7"), 0«(t), and ay (T) is good 
for oil and water droplets in air whenever T > 0.5, we can envisage situations for which 0<k(T) evaluations 
will be adequate. From these numerical examples, we conclude then that whenever /3 < 1 and for r > 0.5, 
the numerical predictions of the zero-gradient and finite-gradient models quantitatively agree to within 10 
per cent. This of course will not be true if /3 > 1. Hence, for the oil and water droplets considered here, the 
predictions of the zero-gradient and finite-gradient models differ significantly only for values of T substantially 
less than 0.25. 

Another quantity which enters the design of combustion chambers is the time required to achieve the 
superheat limit of emulsified fuel droplets; i.e., the time required before the emulsified fuel droplet 
undergoes a "microexplosion." For illustrative purposes, let us assume that the superheat limit T Sh for an 
oil-water emulsification is 260 °C, and the droplet temperature T^q is 20 °C, and the gas temperature T g() is 
1520 °C. The superheat limit that corresponds to the reduced temperature SL = 0. 16. From the numerical 
examples, we find that for pure oil droplets, d oWio = 0.313) = 0.16, 0^(7%, = 0.317) = 0.16, 0,(t? = 
0.233) = 0.16, and av (Tav = 0.286) = 0.16. The corresponding reduced time values at the superheat 
limit for pure water droplets are 

M4 = 1.45) = 0.16, 6 d Mz (= 1.47) = 0.16, 
8 (f = 1.29) = 0.16, and av (T^ = 1.36) = 0.16. 

Using the values of t& = 36.8 ms for oil and ta = 15.6 ms for pure water, we list the superheat times 
predicted by the two models for pure oil droplets; namely, t% = 11.5 ms, tfe — 11.7 ms, fj =8.7 ms and 
ta V = 10.5 ms. The corresponding values for pure water droplets are f§o = 22.5 ms, f& = 22.9 s, tf = 
20.1 ms, and £ av = 21.1 ms. Again, the quantitative differences between the two models for the cases 
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FlGURE 1. The reduced temperature 6 as a function of reduced time T 
predicted by the zero gradient model for oil droplets in air. 

The solid curve with solid dots corresponds to the reduced temperature d oz given by eq (26). The 
dashed curve with open circles corresponds to the reduced temperature &&, given by eq. (31). All of 
the above quantities are dimensionless and the droplet radius is 50 /im. 



Figure 2. The reduced temperature 6 as a function of reduced time r 
predicted by the finite gradient model for oil droplets in air. 

The solid curve with solid dots is the average droplet temperature d av given by eq (36) and the 
dashed curve with open circles is the surface temperature of the droplet 9 t given by eq (34) with TJ = 
1. All of the above quantities are dimensionless and the radius of the droplet is 50 Aim. 




FIGURE 3. The reduced temperature S as a function of reduced time T 
predicted by the zero gradient model for water droplets in air. 

The solid curve with solid dots corresponds to the reduced temperature 0^, eq (26), and the 
dashed curve with open circles corresponds to the reduced temperature Ofo, eq (31). All of the above 
quantities are dimensionless and the radius of the droplet is 50 /Am. 



Figure 4. The reduced temperature d as a function of reduced time T 
predicted by the finite gradient model for water droplets in air. 

The solid curve with dots is the average droplet temperature 6 av , eq (36), and the dashed curve 
with open circles is the surface temperature of the droplet t , eq (34), with tj = 1. All of the above 
quantities are dimensionless and the radius of the droplet is 50 //.m. 
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considered in this paper are not significant. Here the subscripts denote the model and the superscript 
denotes pure oil and the superscript w denotes pure water. 

In summary then, whenever the droplet thermal conductivity exceeds the thermal conductivity of the 
surrounding gas and the reduced time r exceeds 0.5, the quantitative predictions for the reduced 
temperatures of the zero-gradient and finite-gradient models agree to within 10 percent and the superheat 
times agree to within 25 percent for pure oil droplets and to within 15 percent for pure water droplets. 
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7. Appendix A 

In this appendix, we prove that the function A given in eq (29) has no complex roots. 

We first assume that such roots exist and then show that this assumption leads to a contradiction. Let us 
assume that pi and p 2 are nonzero and that A(p$) = where i = 1 or 2. We define the auxiliary functions 

Hi(x) = X" 1 sinhp^ ( A1 ) 

and 

^(X) = X" 1 sihn (Pi) e *p {-<xPi(7l ~ !)}• ( A2 ) 

These auxiliary functions satisfy the differential equations 

p?|a, = 0< X <1 (A3) 






a 2 p 2 Au 2 = 1<X (A4) 

and the conditions that Ui(0) = 0, u 2 ( cc ) = 0, ^i(l) = u 2 {l), and u'i(l) = (3u' 2 (l). The last condition is a 
consequence of the definition of p x . Let f i(x) and v 2 ix) De tne analogues of Ui(x) and u 2 (x) with p x replaced 
by p 2 and consider eqs (A5) to (A8) to be the corresponding versions of eqs (Al) to (A4). For brevity, we do 
not list eqs (A5) to (A8). Equation (A3) is now multiplied by V\, and equation (A7) is multiplied by U\. The 
two are subtracted. The resulting equation is multiplied by \i and then integrated from zero to one to obtain 

(pi - pi) xV>i dx + X 2 ("i"i - uy t ) dx + 2x(u 1 v l - u x v[) dx = 0. (A9) 

Jo 4) Jo 

Performing the same operations on the functions u 2 and v 2 but integrating instead from one to infinity yields 

free /CO 

X 2 u 2 v 2 dx + J X 2 ("2"2 - u 2 v' 2 ) dx+ J 2x(u' 2 v 2 - u 2 v 2 ) dx = 0. (A 10) 

Multiplying eq (A10) by /3, subtracting it from eq (A9), and using the conditions on u u v±, u 2 , and v 2 , we 
obtain the following condition: 

(P? - P?) { | xVi d X + pc? ^ X W 2 d X ) = 0. (All) 
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We now assume that both p l and p 2 are complex and that p 2 = p*. In this case, uf — V\ and u$ = v% so 
that both integrals in eq (All) are of positive definite quantities and therefore necessarily positive. The curly 
brackets then cannot vanish implying that p\ — pf = 0. This criterion cannot be satisfied by pi and p 2 of 
the assumed nature. This yields the contradiction mentioned earlier and is sufficient to rule out the existence 
of complex roots. 



8. Appendix B 

As stated in section 5, the integrals representing the temperature distributions of the droplet and 
surrounding gas, Od(r), r) and 6 g (r) , r) given as eqs (33) and (34), must be evaluated numerically. The same 
remark also holds of course for the average temperature in the droplet, ^(avg, r), given as eq (35). Let us 
denote the denominator of the integrands of each of these three expressions by 

D({) = C?/?? sin 2 € + (f cos f - (1 - 0) sin f) 2 . (Bl) 

When/)(f) becomes small, the integrands become large. 

The values of a 2 /3 2 for the oil and water droplets are from table 2, respectively, 1.13 X 10 -4 and 1.24 X 
10 -5 . Such small values lead to very large integrands whenever f is near one of the minima of D(jj). We may 
approximate the values of ft for which /}(ft) is at one of its minima by the roots of the equation 

ftcosft-(l - #sinft = (B2) 

when a 2 fi 2 is small compared to one. As a consequence of such small a 2 /3 2 , the integrands are highly 
oscillatory. For this reason, Gaussian-type integration formulas are inappropriate. They would require an 
extremely high point formula in order to obtain accurate results. The procedure used consisted of two steps 
which we now outline schematically. The first step is to terminate the upper limit of integration so that 

f /{x)dx = f Ax) dx, (B3) 

Jo Jo 

Where/? is choosen to make the result independent of/? to the desired accuracy. Secondly, a function g(x) 
which behaves \ikef(x) in the neighborhood of its large oscillations is introduced as follows 

f Mdx = I \Ax) - g(*)]dx + f g(*) dx. (B4) 

Jo Jo Jo 

The function g(x) is chosen subject to two restrictions; namely, that /(:*;) — g(x) is a relatively smoothly 
varying function of x, and secondly that the last integral in eq. (B4) may be evaluated analytically. Finally, 

I [/W — g(%)\dx is computed using a finely-spaced Simpson's rule. 

We obtain a function g(x) by summing three expressions derived from the first three roots, f i, %2-> and £3 
of eq (B2). The expression which corresponds to the ith root is the numerator of the particular integrand in 
one of the eqs (33), (34), and (35) evaluated at i; = (j $ with the denominator replaced by 

D($ ^ c?f?& sin & 4- (0 cos ft - ft sin ft) 2 (f - ft) 2 . (B5) 

The initial conditions, 6 d (rj, 0) = and o (r), 0) = 1, were checked using R = 20 For r > 0.25, the 
relative error introduced in using R = 10 instead of R = °° was never larger than 10 -6 for any case. The 
results displayed graphically in section 6 were computed to at least four significant figures. This required 
varying the spacing in the Simpson's rule until the desired accuracy was achieved. The spacing was 
generally of the order of 10 -4 . 
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